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We develop a theory for a generic instability of a Fermi liquid in dimension d > 1 against the 
formation of a Luttinger-liquid-like state. The density of states at the Fermi level is the order param- 
eter for the ensuing quantum phase transition, which is driven by the effective interaction strength. 
A scaling theory in conjunction with an effective field theoy for clean electrons is used to obtain the 
critical behavior of observables. In the Fermi-liquid phase the order-parameter susceptibility, which 
is measurable by tunneling, is predicted to diverge for 1 < d < 3. 
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Landau's Fermi-liquid theory provides a very success- 
ful paradigm in condensed matter physics. By map- 
ping the low-lying excitations in interacting Fermi sys- 
tems onto those of noninteracting ones |]J it explains 
many properties of electrons in solids, including the lin- 
ear temperature (T) dependence of the specific heat, and 
the quadratic T-dependence of the electrical resistivity 
0. In a renormalization-group (RG) context it can be 
understood as the scaling behavior near a stable fixed 
point (FP) that governs the low-T behavior of the system 
[3|. Because of this success, deviations from Fermi-liquid 
(FL) behavior have attracted considerable attention 4]. 
Examples includeparts of the normal phase of high-Tc 
superconductors [HI], heavy-fermion systems [3], and the 
paramagnetic phase of the helimagnet MnSi at low T 7] . 

There are different sources for non-Fermi-liquid (NFL) 
behavior. One is the vicinity of a quantum critical point 
where the low-T behavior is governed by a critical FP 
rather than the stable FL FP, which changes the scaling 
behavior. The electrons couple to the critical soft modes, 
and the NFL behavior is confined to a small region in 
parameter space. This is believed to be the source of 
the observed NFL behavior in heavy-fermion systems . 
Another possibility is the existence of Goldstone modes 
due to a spontaneously broken symmetry and resulting 
long-range order. If the electrons couple to the Goldstone 
modes, NFL behavior can result in an entire phase. This 
has been proposed to explain the behavior of MnSi Q. 
More generic mechanisms for NFL behavior, that do not 
rely on underlying long-range order, are hard to find. It 
is well known that in one-dimensional (1-d) fermion sys- 
tems an arbitrarily small repulsive interaction amplitude 
Kg (we will restrict ourselves to a point-like interaction in 
the spin-singlet particle-hole channel) leads to an insta- 
bility of the FL against a Luttinger liquid (LL) that has a 
vanishing density of states (DOS) at the Fermi level @. 
The LL is characterized by sound-like excitations, and 
its properties, apart from the vanishing DOS, are very 
similar to those of a FL. A natural question is whether 



in dimensions > 1 a similar instability will occur for Kg 
greater than some > 0. Despite substantial efforts, to 
date no description of such an instability has been found. 

There are, however, indications that an instability ex- 
ists. Perturbation theory in the FL phase yields non- 
analytic dependencies on T, or the wave number fc, for, 
e.g., the spin susceptibility and the specific heat coeffi- 
cient IB" 12 1 • For generic d they take the form T"*"^ 
or k'^~^, with multiplicative logarithms in odd d. For 
d = 1, the logarithmic divergencies coincide with the 
well-known perturbative signatures of the LL [l^ . This 
is reminiscent of disordered electrons, where perturba- 
tion theory generically yields a T'^''^^^/^ or fc'^^^ behavior, 
which turns into logT or logfc in d = 2. These pertur- 
bative "weak-localization" effects signalize the instability 
of the disordered FL against an Anderson or Anderson- 
Mott insulator [ij, |l5| . In d = 2 this instability occurs at 
arbitrarily small values of the disorder, whereas in d > 2 
a metal-insulator transition occurs at a nonzero critical 
value of the disorder. It is thus natural to speculate that 
a transition from a FL to a LL can occur in d > 1 . 

In this Letter we construct a theory that describes a 
quantum phase transition from a FL to a NFL state with 
a vanishing DOS at the Fermi level in d > 1. The DOS 
serves as the order parameter (OP) for the transition; 
the FL is the ordered phase, d^ = 1 is the lower criti- 
cal dimension for the transition; fluctuations destroy the 
ordered phase for d < d~ . For d = l-|-e (e^ 1) the the- 
ory is controlled and the critical value of the interaction 
strength is = 0(e^/^). For larger d the critical behav- 
ior is obtained from scaling considerations. Our approach 
is inspired by the nonlinear sigma-model description of 
the classical magnetic Heisenberg transition near d = 2 
[l^ . In many respects our theory is analogous to that 
of the Anderson-Mott metal-insulator transition of dis- 
ordered interacting electrons, for which d~ = 2 as well 
Til 17 1 , even though the transition in disordered systems 



is to a non-standard insulator, while in clean ones, it is 
to a non-standard metal (isl. 
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To identify the DOS as the OP for the FL-to-NFL tran- 
sition we consider a Ward identity that reflects the bro- 
ken symmetry between retarded and advanced degrees of 
freedom in a FL. It relates a two-particle correlation func- 
tion F2 (schematically, (ipip'ipip) , with ip and ip fermion 
fields; a more explicit expression of F2 will be given in Eq. 
p3)) ) to a single-particle function Fi ((V'l/')) and is a gen- 
eralization of a Ward identity first considered for nonin- 
teracting electrons with quenched disorder [l^jliOl- With 
p the center-of-mass wave vector {\p\ ^ k^; we denote the 
Fermi wave number, energy, and velocity by fcp, ep, and 
Up), k the hydrodynamic wave vector (|fc| = k <C fcp), 
and TTic the electron mass it can be written 

(iri„i_„2 + p ■ k/ruc) F2 {p, k; iuj^ , i^nj 

= ^"1 (P, fc;»w„i,za;„2) , (la) 

Fi is proportional to the difference between Green func- 
tions taken at the fermionic Matsubara frequencies iun-i 



and 



For zri„i_„2 



0, fe ^ 0, Fi 



vanishes if w„^w„2 > 0, but is nonzero if < 0. In 

the latter case, and for noninter acting electrons, 

Fi (p, k;iuJni, i^n2 ) oc i sgn (Jl )(5(ep-eF), (lb) 

with tp the single-particle energy-momentum relation. 
For w„j^w„2 < there thus is a family of 4- fermion func- 
tions that diverge in the hydrodynamic limit of vanishing 
ri„j_„2 and k. Taking moments with respect to p yields 
an infinite number of soft modes, provided the density of 
states at the Fermi level, cx (5(ep — ep), is nonzero. 
In the presence of quenched disorder, in contrast, only 
the zeroth moment of Eq. ([Ta|) is soft. These soft modes 
are the Goldstone modes of a spontaneously broken ro- 
tational symmetry in frequency space that is a nonlocal 



generalization of the one considered in Ref. |21 



Equations. ([T]) can be generalized to interacting sys- 
tems along the lines of Ref. [2l| . FL theory ensures that 
the basic structure of the identity is unchanged: The 
symmetry is broken, and Goldstone modes exist, as long 
as the DOS at the Fermi level is nonzero. The nonin- 
teracting DOS, iVp, gets replaced by the physical DOS, 
7V(ep), and the prefactor of the frequency i7„j^_„2 ac- 
quires a FL correction. F2 remains massless, and the 
frequency continues to scale as a wave number. Con- 
versely, a vanishing DOS implies that the symmetry is 
restored and the Goldstone modes have zero weight. If 
this happens, by varying some control parameter, with 
the Fermi level inside the conduction band, then the sys- 
tem will undergo a symmetry-restoring phase transition 
from a FL (ordered phase) to a NFL (disordered phase) 
with the DOS as the OP. In the FL the Goldstone modes 
are all proportional to the basic Goldstone propagator 

V{k, iVL) = A^(ep) ip{iVt/vFk)/k. (2a) 

The explicit form of the function (p depends on the di- 
mensionality. In the limiting case d — ^ 1 one has 

Lp{x) oz\x\/{l + x^). (2b) 



These considerations show that in the ordered phase 
there are soft modes whose frequency scales as a wave 
number, ~ A:: A dynamical exponent z = 1 is asso- 
ciated with the stable FL FP. At a symmetry-restoring 
transition, described by a critical FP, z ^ 1 in general. 

The Goldstone modes are not related to the density 
propagator, which is governed by particle-number con- 
servation. The latter, plus the fact that the thermody- 
namic density susceptibility dn/d)! is expected to show 
no critical behavior (see below), implies that in the den- 
sity propagator one has 17 ~ A;, or z = 1, at both the FL 
and the critical FPs. This is consistent with the fact that 
ri ~ fc at the stable FP that describes a LL in d = 1 . 
Hence there is more than one dynamical exponent: The 
critical dynamical exponent z (2 7^ 1 in general), related 
to the Goldstone modes, and a second dynamical expo- 
nent Zc = 1 related to the charge or density dynamics. 

In what follows, we construct a scaling theory for a 
symmetry-restoring FL-to-NFL quantum phase transi- 
tion where the DOS vanishes. We have also derived an 
effective field theory that allows for an explicit descrip- 
tion of such a transition, the most important aspects of 
which we will sketch at the end of this Letter. 

We start by considering the free energy density /, 
which quite generally satisfies a scaling relation 

f{t,T,h)^b-^'^+'-^ f{tb^/'',TF',hby'^) . (3) 

We have assigned scale dimensions [L] — —1 and [T] — z 
to factors of length and temperature, energy, or inverse 
time (fi= fee = 1), which yields [/] = —d~z for the scale 
dimension of / |22|. b is the RG length rescaling factor. 
h is the field conjugate to the OP, with scale dimension 
[h] = Uh- t is the dimensional distance from the critical 
point, and v = l/[t] is the correlation length exponent. 
For the OP density N = -{df/dh)/T Eq. ^ implies 



N{t,T) = b-'^+y^ N{tb^^'',Tb') 
= b-'^+'- N{tb^'\T¥') . 



(4a) 
(4b) 

Equation (j4b|) reflects the fact that N is the DOS, which 
scales as an inverse energy times an inverse volume. 
Hence yu = z. At T = and at criticality, respectively, 
the OP vanishes as a power law, 

N{t, T^O)(xt^ , N{t^ 0, T) (X T^'^-'-)/'- . (5a) 

with a critical exponent 

I3 = iy{d~z). (5b) 

The corresponding results for the specific heat coefficient 
7 are obtained from Cy ^ = -Td^f/dT^. The 
scaling behavior of 7 is the same as that of the DOS: 



7(t,T) = 6-'*+^7(i6i/'',T6"), 



(6) 



We next consider the OP susceptibility x = dN/dh as a 
function of t, T, and the wave number k. In general, 

x{t,T;k)^b^-^x{tb'/'',Tb^,kb), (7a) 
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which defines the exponent rj. At T = and at criticality, 
respectively, the homogeneous OP susceptibihty diverges: 

X(i,r = 0,fc = 0) cx t-^ , 7 = ^.(2-7?), 

x(t = o,r,fc = o) cx r-(2-'))/^ (7b) 

From Eqs. ([3]), (HI, and ([7a|) we find the exponent relation 

z = {d-T] + 2)/2. (8) 

This implies that there are only two independent criti- 
cal exponents, e.g., v and z (see, however, the remark 
above regarding multiple exponents z) rather than three 
as is generally the case at a quantum critical point 
(iij . For dn/d^ we expect no critical behavior since it 
does not show the perturbative nonanalyticities that are 
precursors of the critical behavior of other observables 
jiol [23 |. The scaling behavior of the electrical conduc- 
tivity a = Dcdn/d^ is therefore given by that of the 
charge diffusion coefficient Dc, which scales as a length 
squared divided by a time. Since describes the charge 
or density dynamics the relevant dynamical exponent in 
this context is Zc = 1- We thus have [25| 



(9a) 



If z < 1 (see below) this yields for the electrical resistivity 
p — 1/a at criticality 



p{t = 0,r) cx T. 



(9b) 



The preceding scaling predictions all pertain to the 
critical FP. Also of interest are the OP and the OP sus- 
ceptibility in the ordered phase, \t\ — oo, where rj — d, 
which implies z — 1. From Eq. (j4bp we have 



N{\t\ = c»,T) cx const. + r' 



d-l 



(10) 



This is one example of the perturbative nonanalyticities 
mentioned above. The same power law holds at T = 
as a function of the distance lu from the Fermi surface: 
N{T — 0,U!) cx const. + ui'^^^. It is analogous to the 
Coulomb anomaly in disordered systems, where N{T — 
0, uj) cx const. -I- w''*"^-'/^ [2^. The latter is a precursor of 
the quantum phase transition in disordered systems (the 
Anderson- Mott transition [3l), where the DOS vanishes 
and serves as an OP The current theory suggests 

that an analogous statement holds in clean ones. For the 
OP susceptibility we find from Eq. ([7a| 



X(K| =cx),T,fc) 



k'^-^f^iT/k) cx r/fc3-^. 



(11) 



In the second relation we have used the result of an 
explicit calculation [2I], which yields f^{x — >■ 0) cx x. 
This divergence of the OP susceptibility, or the 2-point 
local-DOS correlation function, which is observable by 
tunneling experiments, is a consequence of the Gold- 
stone modes. It is analogous to the l/fc'*"'' divergence 



of the longitudinal susceptibility in the ordered phase of 
a Heisenberg ferromagnet [2^. For a 2-d FL it predicts 
a 1/fc divergence with a prefactor that is linear in T. 

The preceding scaling considerations are expected to 
be valid between the lower critical dimension d~ — I 
and some upper critical dimension d^ . Equation (jlip 
suggests (i+ = 3, but this requires further corroboration. 
For d > d^ one expects the critical behavior to be mean- 
field like and governed by a Gaussian FP. 

We now sketch the derivation of an effective field the- 
ory that allows for an explicit description of a quantum 
phase transition of the type we have discussed above. A 
complete account will be given elsewhere [2^. This ef- 
fective theory is in the spirit of the matrix field theories 
that were pioneered by Wegner 2^ S^, and generalized 
by others Sll, for disordered systems. We consider a 
fermionic action and define electron bispinors 

■r]n{x) = {4>n1:{x),'ipnl{x),^/jnl{x),-1pnt{x)) / V2 (12) 

where ip and ip are fermionic fields with Matsubara fre- 
quency index N and spin projection fi, as well as adjoints 
rinix) = Cr]n{x) with C = ici (8)cr2, where ai,2 are Pauh 
matrices. We confine the tensor product ri^{x)®rjrn{y) to 
a spin-quaternion- valued bosonic matrix field Qnmix,y) 
by means of a Lagrange multiplier A„m(a;, y). The Ward 
identity then takes the form of Eqs. ([1]) with 



xtr( 



^ (trg„,„,(p + fe/2,p- fc/2) 
>nin.(p-fe/2,P + fc/2)) (13) 



where a;„jCJ„2 < 0. This identifies qnm{Pi,P2) = 
Q{—nm)Qnm{PiiP2) as the Goldstone modes. The cor- 
responding elements A of the Lagrange multiplier matrix 
field A are also soft modes. The electron-electron inter- 
action couples q to Pnm{Pl,P2) = Q{nm) Q„rn{Pl , P2) , 

and integrating out P, and the corresponding part of A, 
generates terms to all orders in q and A. An analogous 
procedure in the presence of quenched disorder provides a 
perturbative derivation, order by order in powers of q, of 
the generalized nonlinear sigma-model for the Anderson- 
Mott transition problem [l3l ■ We have derived the 
complete action to order g*, which suffices for a 1-loop 
RG calculation. The effects of A can be absorbed into 
diagram rules. The same method can be used to system- 
atically derive higher-order terms. 

This effective theory can be analyzed by RG methods 
in. d = 1 -I- e in analogy to the disordered case in d = 2 -I- e 
[l5| . The 1-point function is proportional to the DOS: 



P'^^^ = (trQ„„(a;,a;))| 



cx Niep) = NpZ^^^. (14) 



Physically, Z^/^ = (1 + SZy^^ is the physical DOS nor- 
malized by the bare or free-electron DOS; technically, it 
is the field-renormalization constant. It is related to, but 
not the same as, the residue of the pole in the Green 
function. At 1-loop order one finds that SZ is negative. 
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logarithmically divergent in d ~ 1, and proportional to 
1/e in d = 1 + e. As a function of Kg it is of 0{K^) for 
small Ks- In a naive extrapolation the DOS thus vanishes 
at a critical value = 0(e^/^). The 2-point function 

^'^^ = {qnin2{kl, k2) Qnsniiks, ki)) (15a) 

has a contribution proportional to (5„i„3 (5„2„4 that con- 
stitutes the Goldstone propagator V, Eqs. In c? = 1 

V{k,in) = Z H\n\/{k^/G^ + H^n^). (ISb) 

The bare values of G (which is the loop expansion pa- 
rameter) and H are l/vpNp and Np, respectively. An 
explicit calculation finds no singular renormalizations of 
G and H at 1-loop order. Structural considerations con- 
firm this and show that at 2-loop order there is a singu- 
lar renormalization of G due to insertion diagrams; an 
inspection of skeleton diagrams will require a 2-loop cal- 
culation. This strongly suggest a critical FP at 2-loop 
order with a FP value of the renormalized coupling con- 
stant g — b^G given by g* — 0(e^/^). Choosing the 
independent exponents to be v and z this leads to 

i^ = l/2e + 0(l) , z = l + 0(e). (16) 

The 0(e) term in z requires a 2-loop calculation; an ed- 
ucated guess is as follows. In the bare theory, G^H cx 
rrLc/n, with n the electron density. This quantity one 
does not expect to be renormalized, so H ^ G^^ ^ b^^'^, 
or z = 1 — e [sij. Hence z < Zc, which implies Eq. (j9b[) . 

We stress that the RG theory, although it relies on 
an e-expansion about d = 1, is not tied to the special 
properties of 1-d fermion systems [9]. It only requires 
Goldstone modes with fc ^ f2, which exist in any d > 1, 
and we can apply a RG scheme that works in a fixed 
dimension [s^l d > 1. 
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